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Abstract 

A search for narrow-width resonances that decay into electron+jet or neu- 
trino+jet has been performed with the ZEUS detector at HERA operating at 
center-of-mass energies of 300 and 318 GeV. An integrated e + p luminosity of 
114.8 pb" 1 and e~p luminosity of 16.7 pb _1 were used. No evidence for any 
resonance was found. Limits were derived on the Yukawa coupling, A, as a func- 
tion of the mass of a hypothetical resonance that has arbitrary decay branching 
ratios into eq or vq. These limits also apply to squarks predicted by i?-parity- 
violating supersymmetry. Limits for the production of leptoquarks described by 
the Buchmuller-Rtickl-Wyler model were also derived for masses up to 400 GeV. 
For A = 0.1, leptoquark masses up to 290 GeV are excluded. 
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1 Introduction 



Many extensions of the Standard Model (SM) predict the existence of particles carrying 
both baryon and lepton numbers, such as leptoquarks (LQs) |I] or squarks (in i?-parity- 
violating supersymmetry) j2]. In ep collisions at HERA, such states may be produced 
directly through electron 1 -quark fusion, with subsequent decay into electron and quark 
or neutrino and quark, yielding peaks in the spectra of the final-state lepton-jet invariant 
mass, Mij. This paper presents a search for such resonant states. 

The only significant backgrounds to high-mass resonance production arise from neutral 
current (NC) and charged current (CC) deep inelastic scattering (DIS), as illustrated in 
Figs, and b. While a resonance with mass below the HERA center-of-mass energy, 
y/s, would give rise to narrow peak in the My spectrum, the backgrounds from NC and 
CC fall rapidly at high mass due to the dependence of the cross section on Q 2 , the 
virtuality of the exchanged bosons, and to the sharply falling valence-quark density at 
large Bjorken x. The variable 8*, the lepton scattering angle in the lepton-jet center- 
of-mass frame, can be used to reduce the DIS backgrounds. The decay of a resonance 
results in an angular distribution different from those produced by SM processes: a scalar 
resonance, for example, will have a flat distribution in cos#*, while NC DIS events follow 
approximately a 1/(1 — cos#*) 2 distribution. 

The data were used to investigate the production of resonances. In the absence of a narrow 
resonance signal, limits can be set on the production of resonances with masses below y/s 
and with the assumption of narrow width, as shown in Fig. [TJ;. In addition, specific limits 
on the production of Buchmuller-Riickl-Wyler (BRW) LQs pQ with masses both below 
and above y/s can be obtained. Because backgrounds fall sharply at large My, the data 
are sensitive to LQs via exchange terms, shown in Fig. ^i, as well as interference terms 
with SM processes. The narrow-width assumption is not necessary for setting limits in 
the BRW model. 

The present results supersede previous analyses [HE]. They are based on all data collected 
by the ZEUS experiment in the period from 1994 to 2000. After 1998, the HERA center- 
of-mass energy was increased from 300 to 318 GeV. 

2 Experimental conditions 

The ZEUS detector is described in detail elsewhere [3]. The main components used in 
the present analysis are the central tracking detector (CTD) and the uranium-scintillator 

1 Unless otherwise specified, 'electron' refers to both positron and electron and 'neutrino' refers to both 
neutrino and antineutrino. 
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sampling calorimeter (CAL). 

Charged particles are tracked in the CTD j^J, which operates in a magnetic field of 
1.43 T provided by a thin superconducting solenoid. The CTD consists of 72 cylindrical 
drift chamber layers, organized in nine superlayers covering the polar-angle region 15° < 
9 < 164°. The transverse-momentum resolution for full-length tracks is ct(pt)/Pt = 
0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The CAL [Zj consists of three parts: the forward 2 (FCAL), the barrel (BCAL) and the 
rear (RCAL) calorimeters. Each part is divided into modules, which are subdivided trans- 
versely into towers and longitudinally into one electromagnetic section (EMC) and either 
one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). The smallest sub- 
division of the calorimeter is called a cell. The CAL energy resolutions, as measured under 
test-beam conditions, are a(E)/E = 0.18 /y/E for electrons and a(E)/E = 0.35/ y/E for 
hadrons, with E in GeV. The timing resolution of the CAL is better than 1 ns for energy 
deposits greater than 4.5 GeV. Further performance parameters of the CAL relevant for 
this study have been discussed in previous publications OH]- 

The forward plug calorimeter (FPC), a lead-scintillator sandwich calorimeter, was in- 
stalled in the 20 x 20 cm 2 beam hole of the forward CAL (FCAL) of the ZEUS detector 
in 1998. Although the FPC information was not used in this analysis, the impact of its 
material on the CAL response to forward jets was extensively studied. 

The luminosity, which was measured [S] from the rate of the bremsstrahlung process 
ep — > epj, has an uncertainty of 1.6% to 2.25%, depending on the running periods. All 
ZEUS data collected from 1994 to 2000 are listed in Table ^ an d were used for this 
analysis. 

3 Monte Carlo simulation 

Production and decay of resonances were simulated using PYTHIA 6.1 (HI CD], which 
takes into account the finite width of the resonant state, but includes only the s-channel 
diagrams. Initial- and final-state QCD radiation from the quark and the effect of LQ 
hadronization ^2] before decay as well as the initial-state QED radiation from the electron 
are also taken into account. 

Standard Model NC and CC DIS events were simulated using the HERACLES 4.6.2 [H] 
program with the DJANGO 6 version 2.4 [TJ] and DJANGOH 1.1 [T3] interfaces to 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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the hadronization programs. Radiative corrections for initial- and final-state electroweak 
radiations, vertex and propagator corrections, and two-boson exchange were included. 
The hadronic final state was simulated using the MEPS model in LEPTO 6.5 [T3], which 
includes 0(as) matrix elements and higher-order QCD radiation. The CTEQ5D parton 
distribution function (PDF) [T^j was used in evaluating the SM cross sections. 

The largest uncertainty in the NC and CC cross sections is due to the uncertainties in 
the parton densities of the proton. The PDFs at high Bjorken x are determined primarily 
from measurements made in fixed-target DIS experiments. At x — 0.6, corresponding to a 
lepton-jet mass of 230 GeV, the cross-section uncertainty due to the PDF uncertainty fH] 
is ~ 6% for NC and ~ 4 (10)% for e~u (e + d) CC reactions, where u and d refer to the u 
and d valence quarks, respectively. 

The generated events were passed through the GEANT 3.13-based [T7] ZEUS detector- 
and trigger-simulation programs j3] . They were reconstructed and analysed by the same 
program chain as the data. 

4 Resonance search 

Events from a hypothetical resonance decaying into eq (uq) have a topology identical to 
DIS NC (CC) events. Events originating from high-mass resonances are expected to have 
high transverse energy, at least one jet, and either an identified final-state electron or 
large missing transverse energy. The lepton-jet invariant mass was calculated as 



where Ei is the energy of the outgoing lepton, Ej is the energy of the jet, and £ is the 
angle between the lepton and jet. In final states containing multiple jets, the jet with the 



4.1 ep —> eX topology 
4.1.1 Event selection 

Events with the topology ep — > eX, where X contains one or more jets, were selected 
using the following criteria: 

• the Z coordinate of the reconstructed event vertex was required to be in the range 
\Z\ < 50 cm, consistent with an ep collision; 




largest transverse momentum, P^, was used. 
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• the total transverse energy, E?, was required to be at least 60 GeV. This removes the 
bulk of the SM NC background; 

• an identified electron ^H] was required with energy E' e > 25 GeV located either in 
FCAL or BCAL, corresponding to an electron polar angle, 8 e < 126°. Electrons 
impacting the BCAL within 1.5 cm of a module edge, as well as electrons impact- 
ing between the FCAL and the BCAL, as determined by tracking information, were 
discarded to ensure that the resolutions were well understood; 

• at least one hadronic jet with transverse momentum Pj, > 15 GeV, obtained using the 
longitudinally invariant cluster algorithm JH] i n inclusive mode [20], was required. 
The centroid at the FCAL face of the highest-Py jet was required to be outside a box of 
60 x 60 cm 2 centered on the proton beam, in order to ensure good energy containment 
and to reduce the systematic uncertainties due to the proton remnant. 

The acceptance, mass shifts, and resolutions for resonant lepton-quark states were calcu- 
lated from the LQ MC. After these cuts, the acceptance for scalar resonances was ~60%, 
depending weakly on the mass. The acceptance for vector resonances was ~60% below 
200 GeV, decreasing to ~40% at 290 GeV. These differences in the acceptances for 
scalar and vector resonances are due to the different decay angular distributions. The 
mass resolution, determined from a Gaussian fit to the peak of the reconstructed mass 
spectrum, fell from 6% to 4% as the resonant mass increased from 150 to 290 GeV. The 
peak position was typically lower than the generated mass by 1%. The resolution in cos#* 
near | cos#*| = 1 had a Gaussian width of 0.01, degrading to 0.03 with decreasing | cos0*|. 

4.1.2 Search results 

After the above selection, 21 509 events were found, compared to 21445 ± 1 288 expected 
from the NC MC simulation and the evaluation of its systematic uncertainties (see Sec- 
tion ^T3J). The measured distributions of the total transverse energy, Et, are compared 
to the simulation in Figs. |2K and|2t, where the e + p and e~p samples are shown separately. 
Also shown are E — Pz (Figs. I2t>,f), where the E and Pz are summed over the final-state 
electron and all the hadrons, electron transverse momentum (P<f) (Figs. |2t,g) and jet 
transverse momentum {Pj) (Figs. |2Ji,h). Good agreement is seen between the data and 
the SM NC simulation for all of these spectra. 

Figures El and |U show the M e j spectra for e + p and e~p data, respectively. The upper plots 
show the spectra with and without the cut cos#* < 0.4, while the lower plots show the 
ratio of the observed spectrum to SM expectations with no cos#* cut. The data are well 
described by the NC MC. 
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An excess of data events relative to SM expectations was seen in data from the 1994-97 
period [3]. For M e j > 210 GeV, 24.7 events were expected and 49 were observed. No 
such excess is seen in the more recent data. The increase in the proton beam energy 
from 820 to 920 GeV in 1998 had no significant effect on the mass reconstruction or 
the signal acceptance. The effect of the addition of the FPC was studied extensively; 
uncertainties in its simulation were found not to affect the conclusions of the present 
analysis. Combining all the e + p data, 104 events were observed with M e j > 210 GeV, in 
good agreement with the SM expectation of 90 ±15. The systematic uncertainty on the 
expectation is discussed in the next section. The excess seen in the earlier data sample 
must be ascribed to a statistical fluctuation. 



4.1.3 Systematic uncertainties 

The uncertainty on the expected number of events from NC DIS process was investigated. 
The dominant sources were uncertainties in: 

• calorimeter energy scale, of 1% for BCAL electrons, 2% for FCAL electrons and 2% 
for hadrons. This led to an uncertainty of 4 (12)% in the NC expectation for M e j=100 
(220) GeV; 

• the simulation of the hadronic energy flow including simulation of the proton remnant 
and the energy flow between the struck quark and proton remnant. Tests based on 
the SM MC samples yielded variations of the NC background of about 10% for masses 
above 220 GeV; 

• the energy reponse of FCAL towers closest to the beam from the differences observed 
between data and simulation. Tests based on the SM MC samples showed variations 
of the NC background of less than 5% for M ej =220 GeV; 

• the parton densities, as estimated by Botje [IB], which gave an uncertainty of 5% for 
M ej =220 GeV. 

Other uncertainties were investigated and found to be small compared to the above items. 
They are the simulation of the electron-energy resolution, the electron-finding efficiency, 
the jet-position reconstruction, the luminosity determination and the simulation of the 
vertex distribution. The overall systematic uncertainties on the background expectations 
result from summing the contributions from all these sources in quadrature and are shown 
in Figs. an d as the shaded bands. 
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4.1 .4 Significance analysis 



To quantify the level of agreement of the M e j spectra between the SM MC and the data, a 
significance analysis was performed using a sliding mass window of width 3a(M e j), where 
a(M e j) is the mass resolution discussed in Section l4.1.1l The number of events, /i, expected 
from the SM and the number of observed events, N, were compared in each window, and 
the probability of observing N or more events while expecting /i was calculated as 



In the e + p data sample, a minimum probability, P m i n , of 2.6xl0~ 3 was found at mass 
121 GeV, where 2 575 events were observed while 2 436 were expected inside the sliding 
window. A large number of MC experiments were then performed, taking into account 
the systematic uncertainties on the SM expectations, and the P m in distribution was de- 
termined. In 35% of all simulated experiments, the value of P m in obtained was less than 
that found in the data. With the additional cut cos 8* < 0.4 applied to the sample, a 
P m i n of 1.9xl0 -2 was found in the e + p data sample at a mass of 158 GeV. A value of 
Pmin less than that found in the data was observed in 41% of the generated experiments. 
The same test was also done on the e~p data samples. The results are summarized in 
Table 121 These observations show that the data are compatible with the SM expectation, 
and there is no evidence for a narrow resonance in the eq channel. 

4.2 ep — ► vX topology 
4.2.1 Event selection 

Events with the topology ep — > uX, where X contains one or more jets, were selected by 
the following cuts, similar to those used in the CC cross-section measurement [2*T] : 

• the Z coordinate of the reconstructed event vertex was required to be in the range 
\Z\ < 50 cm. The event vertex was reconstructed either using the tracks measured in 
the CTD (for events with large 70 [21], the hadronic scattering angle of the system 
relative to the nominal interaction point) or from the arrival time of the particles 
entering the FCAL (for events with small 70, i.e. outside the CTD acceptance); 

• the missing transverse momentum, Pt > 20 GeV, as measured in the calorimeter; 

• y < 0.9, where y was calculated from the longitudinal momentum measured in the 
calorimeter: y = (E — Pz)/2E e , where E e =27.b GeV is the electron beam energy. This 
cut discards events in which the kinematic variables were poorly reconstructed; 




(1) 



k=N 
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• NC events were removed by discarding events with identified electrons; 

• at least one jet was required with P T > 10 GeV, where jets were reconstructed as in 
the ep — > eX topology, and the centroid at the FCAL face of the highest-Pr jet was 
required to be outside a box of 60 x 60 cm 2 centered on the proton beam. 

The neutrino energy and angle were calculated by assuming that Pt and missing E — 
Pz were carried away by a single neutrino. Monte Carlo simulations of resonant-state 
production indicated that the neutrino energy, E u , and polar angle, 9 U , were measured 
with average resolutions of 16% and 11%, respectively. The average systematic shift in 
E u was less than 2%, while the shift in 9 V was less than 1%. After all these cuts, the 
acceptance was ~ 60%, depending weakly on the mass. 

The invariant mass of the z^-jet system, M U j, was reconstructed using the sum of the 
4-momenta of the neutrino and the highest-Pr jet in the event. The shift and resolution 
of the invariant mass were studied by reconstructing it for the LQ MC events and fitting 
the mass peak with a Gaussian function. The resulting mass shift was below 1% for LQ 
masses between 150 and 290 GeV, with the resolution varying from 8% to 7%. 

4.2.2 Search results 

After the selection, 2 536 events were found, compared to 2 587 ± 217 expected from the 
CC MC simulation and the evaluation of its systematic uncertainties (see Section T4 . 2. 3|) . 
The distributions of Pt are compared for data and simulation in Figs. and |U, where 
the e + p and e~p samples are shown separately. Also shown are E — Pz (Figs. Efcsf), 
where the E and Pz are summed over the final-state hadrons, E v (Figs. Et,g) and P T 
(Figs. 01, h). Reasonable agreement is seen between the data and the SM CC simulation 
for all of these spectra. A small excess compared to MC is observed for the e + p data at 
large Pt (also reflected in the P T and E v distributions). 

Figures El and [7| show the M v j spectra for e + p and e~p data, respectively. The upper 
plots show the spectra with and without a cut of cos 8* < 0.4, while the lower plots show 
the ratio of the observed spectra to SM expectations with no cos#* cut. The data are 
reasonably well described by the CC MC. 

4.2.3 Systematic uncertainties 

The uncertainty on the predicted background from SM CC DIS processes was investigated. 
The dominant sources of uncertainty, which are similar to those described in Section l4.1.3l 
for the ep — > eX case, arise from uncertainties in: 
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• the hadronic energy scale, of 2%, which leads to an uncertainty of 2 (10)% for M V j=100 
(220) GeV; 

• the simulation of the energy deposited in the FCAL regions closest to the forward 
beam pipe, which leads to an uncertainty of ~ 7% for M^=220 GeV; 

• the parton densities, as estimated by Botje [TH], giving 9% and 4% uncertainties for e + p 
and e~p, respectively, for M U j =220 GeV. The PDF uncertainties were also estimated 
using the ZEUS NLO fit [22] with similar results. 

Other uncertainties include jet-position reconstruction and luminosity determination, 
which were small compared to those given above. The overall systematic uncertainties on 
the background expectations were obtained by summing the contributions from all these 
sources in quadrature, and are shown in Figs. Eb and Eb as the shaded bands. 

4.2.4 Significance analysis 

The same significance analysis as applied to the M e j spectra was performed on the M v j 
spectra. The results are summarized in Table H3 The observations are again compatible 
with the SM expectation, and there is no evidence for a narrow resonance in the vq 
channel. 

5 Limits on the production of resonant states 

Since no evidence was found for a narrow resonance in the lepton-jet mass spectra, limits 
were set on the Yukawa coupling A and the production cross section of the eight types of 
resonant states listed in Table |H with masses below \fs. The resonance decay width was 
assumed to be small, and the s-channel production was assumed to be dominant, so that 
the resonance exchange and interference contributions were neglected. 

The limits were set using a likelihood technique involving the observables My s and cos 8*. 
The variable My s is the invariant mass of the lepton-jets system and was calculated from: 

M ljs = ^2E e (E + P z ) ljs , 

where E e = 27.5 GeV is the electron beam energy, (E + Pz)u s was determined using the 
lepton and all jets in the event satisfying P£ > 15 GeV (P^ > 10 GeV in ep —>■ vX channel) 
and 7fc < 3 0111 . 

The Mij S method rather than the My method was used in the limit setting because it 
has better resolution on the reconstructed resonance mass for LQ — > eq events. Using 
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the PYTHIA MC described in Section it was found that the resolution on My s for the 
LQ — > eq events varied from 3% to 2% as the mass varied from 150 to 290 GeV, while 
that on Mij varied from 6% to 4%. For the LQ — > vq events, My s and My have similar 
resolutions. On the other hand, the My s method assumes that the final state consists 
only of the LQ decay products and the proton remnant. It is also affected by the QED 
initial-state radiation (ISR) of a photon from the incoming electron. Therefore, the My 
method is more general and better suited for a resonance search. This is discussed in 
more detail elsewhere [Hill]. 

In the limit-setting procedure, the width was assumed to be small compared to the de- 
tector resolution, so that the narrow-width approximation (NWA) is valid. In the NWA, 
the total cross section of the production of a single state is described at Born level pQ by 

^ = (J+l)lAMio,<), (2) 

where q(xo, L) is the initial-state quark (or antiquark) momentum density in the proton 
evaluated at xq = M|_/ s and virtuality scale M 2 q , and J is the spin of the state. 

The effect of QED ISR on the resonance production cross section was evaluated. It varies 
from —4% at resonance mass 100 GeV to —25% at mass 300 GeV and depends weakly on 
the resonance type. This was taken into account when setting the limits on A. The next- 
to-leading-order (NLO) QCD corrections, the so-called .fT-factors, including the vertex 
loop corrections, gluon radiation from the leptoquark or the quark and other higher-order 
diagrams, have been evaluated for the scalar resonances [2*3Tl24| . The iT-factor varied from 
1.17 (1.17) to 1.15 (1.35) for an eu (ed) resonance with mass varying from 100 to 300 GeV. 
However, no calculation is available for the vector resonances, so for consistency, no NLO 
QCD correction was applied to either scalar or vector resonance. 

For each of the running periods listed in Tabled the Mij S -cos6* plane for 150 < My s < 
320 GeV was divided into bins, labelled i, for each of the eq and vq data samples. For 
resonant states with vq decays, both eq — > eX and eq — > vX samples were used, while for 
resonant states decaying only to eq, only eq — > eX samples were used. The upper limit 
on the coupling strength, Ai imit , as a function of M eq , was obtained by solving 3 

A limit OO 

J dX 2 L(M eq , A) = 0.95 J d\ 2 L(M eq ,X), 
o o 

where L is the product of Poisson probabilities of all cos#*-My s bins convoluted with 
3 With the Bayesian prior assumption of a uniform A 2 distribution. 
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Gaussian distributions for the main systematic uncertainties. The L was calculated as 



where j denotes the source of systematic uncertainty and 5j corresponds to the variation 
of the j th systematic parameter in units of the nominal values quoted in Sections 14.1.31 
and 14.2.31 The index % labels the bin in cos 8*-Mij S and iVj gives the number of events 
observed in that bin. The variable /4, which denotes the expected number of events in 
bin i after the effect of the systematic variations, was calculated as 



where /x^, which depends on the M eq and A, is the number of expected events in bin i with 
no systematic variation and gives the fractional variation of /ij under the nominal shift 
in the j th systematic parameter. This definition of reduces to a linear dependence of 
on each Sj when 8j is small while avoiding the possibility of /4 becoming negative which 
would arise if was defined as a linear function of the 5j's. 

Figure |H1 shows the A limits as a function of f3 eq and j3 vq , the branching ratios for LQ — > eq 
and LQ — > vq, respectively, and as a function of the resonance mass. These limits were 
obtained for the four scalar resonant states listed in Table |U and for A = 0.1 and A = 0.3, 
where the latter (A ps ^Ana>Ew) corresponds to the electroweak (EW) coupling. The 
equivalent plots for vector resonant states are shown in Fig. EH For the e + u and e~d 
resonances vq decays are forbidden by charge conservation. The e~u and e + d resonances 
can decay to either eq or vq. The eq+vq limits, which assume f3 vq + (3 eq = 1, are largely 
independent of the branching ratio, and are typically in excess of 285 GeV for the EW 
coupling. The limits obtained using only the eq (or the vq) data set allow for decay modes 
other than eq and vq, so the eq and vq limits are applicable to a wider range of physics 
models than the combined eq+vq results. Similar limits have been presented by the HI 
collaboration for scalar e~d and e~u resonances |25| . 

For comparison, the limits on scalar resonances reported by the D0 experiment [2E] at 
the Tevatron are shown in Fig. |H| For a scalar resonant state with f3 eq =l, the D0 and 
CDF [22] limits are 225 GeV and 213 GeV, respectively, leading to a combined Tevatron 
limit of 242 GeV [2H] • These limits are independent of both coupling and quark flavor, 
but degrade (in the case of D0 for example) to 204 GeV (98 GeV) as [3 eq decreases from 
100% to 50% (0%). 

The limits presented in Fig. |H1 apply to squarks with ftp couplings to eq. For example, the 
A-limit contours on S e +d with the decay S e +d — > e + d (the dashed curves in Fig. [Sfc) apply 




/4 = ^n( 1 + a ^' 
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to Uj squarks with coupling Ayi and subsequent decay Uj — > e + d, where the subscript 
j indicates the squark generation. The limit curves on S e - U with decay S e - U — > e~u 
(the dashed curves in Fig. |H^) and with S e - U — > vd (the dotted curves in the same plot) 
apply to dj squarks with coupling Anj and subsequent decays dj — > e~u and dj — > vd, 
respectively. With the e~u and vd channels combined and with (3 e - u = f3 u d = 0.5/3, the 
limit for squark dj on \nj\/]3 is 0.1 for mass 276 GeV and 0.3 for mass 295 GeV. In 
this case, the i? p -conserving decay branching ratio of dj — > XkQj is 1-/?, where Xk is the 
gaugino. 

The limit on the resonant-state production cross section, cri imit , was calculated by using 
the NWA as shown in Eq. (j2J) with A = Ai; m i t . Figure El shows the results for the scalar 
and vector resonant states. Limits on the calculated production cross section with only eq 
data samples (assuming (3 eq = 0.5), with only vq data samples (assuming (3 uq = 0.5) and 
with all the data samples combined (assuming (3 eq = (3 vq = 0.5) are shown as a function 
of M eq . Usually the limit becomes stricter after combining the eq and vq data samples. 

6 Limits on BRW leptoquark model 

The two-dimensional (My s -cos 8*) likelihood method described above was also used to 
set limits on the Yukawa coupling, A, of the BRW LQs. The full LQ cross section was 
used, including the s-channel and w-channel contributions and interference with DIS. The 
difference of the limits thus obtained to those obtained by using NWA is negligible below 
the LQ mass, M LQ , of 280 GeV. 

The coupling limits of the 14 BRW LQs listed in Table 03 were calculated as a function 
of Mlq. The presence of leptoquarks with Mlq > y/s would affect the observed mass 
spectra, particularly at high mass, mainly due to n-channel exchange and the interference 
effects. Therefore, it is possible to extend the limits beyond \/s. Coupling limits were 
calculated for masses up to 400 GeV. 

The combined search in the LQ — > eq and LQ — > vq topologies, when applicable, produces 
more stringent limits on leptoquark production. An example is shown in Fig. ITTk for the 
LQ produced through e + d fusion, which decays to e + d and vu with equal probability, 
and for the Sq LQ produced through e~u fusion, which decays to e~u and vd with equal 
probability. The limits obtained with only the eq data samples, only the vq data samples 
and the combined samples are shown for each LQ. If a coupling strength A = 0.3 is 
assumed, the production of an Vq LQ is excluded up to a mass of 284 (286) GeV using eq 
{vq) samples only, while it is excluded up to 386 GeV with the combined samples. For the 
Sq, the mass limit is 298 (301) GeV using eq (vq) samples only and the combined mass 
limit is 351 GeV. 



11 



Also shown in Fig. Eh are the limits on the Sfr LQ from £>0, OPAL [H] and L3 [201 • The 

most stringent limits in the high-mass region (above 200 GeV) from the LEP experiments 
come from an indirect search through the t-channel LQ exchange between the incoming 
electron and positron. 

Figure IT21 shows the coupling limits on the scalar and vector BRW LQs with F = and 
F = 2, respectively, where F = 3B + L is the fermion number of the LQ and B and L are 
the baryon and lepton numbers, respectively. The limits on A^q at M^q = 400 GeV range 
from 0.3 to 1.0. In general, present results are significantly better than LEP limits below 
300 GeV and comparable above 300 GeV. The HI collaboration has presented similar 
limits on F = |H5 and F = 2 |2Sj LQs. The excluded mass regions for BRW LQs with 
A = 0.1 and with A = 0.3 are summarized in Table |3 They range from 248 to 290 GeV 
for A = 0.1 and from 273 to 386 GeV for A = 0.3. 

7 Conclusions 

The total ep data recorded by the ZEUS experiment at HERA were used to search for the 
presence of a narrow-width resonance decaying into a lepton and a jet, with the final-state 
lepton being either an electron or a neutrino. The data samples include 16.7 pb _1 of e~p 
and 114.8 pb" 1 of e + p collisions. No evidence was found for resonance production, either 
in the eq or vq topology. Limits were set on the coupling strength of resonant states 
that could decay in these topologies using a two-dimensional likelihood analysis. With 
the combined eq and vq topologies, a scalar e~u (e + d) resonant state is excluded up to 
a mass of 275 (265) GeV for coupling strength A = 0.1. The combined limits depend 
very weakly on the resonance-decay branching ratio. Limits on the coupling strength of 
Buchmiiller-Riickl-Wyler-type leptoquarks with masses up to 400 GeV are also presented. 
The excluded mass regions depend on BRW LQ type and the coupling. At A = 0.1, they 
range from 248 to 290 GeV, which are the most stringent limits available to date. 
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Period 


E p (GeV) 


E e (GeV) 


(GeV) 


e charge 


luminosity ( pb 1 ) 


3lumi 


94-97 


820 


27.5 


300 


e+ 


48.5 


±1.6% 


98-99 


920 


27.5 


318 


e~ 


16.7 


±1.8% 


99-00 


920 


27.5 


318 


e+ 


66.3 


±2.25% 



Table 1: The characteristics of the three data samples used for the present study; 
o~iumi is the uncertainty of the measured luminosity. 





No cos 9* cut 


cos6* <0.4 




e + p 


e p 


e + p 


e p 


p . 

1 mm 


0.26% 


0.032% 


1.9% 


4.3% 


M ej (GeV) 


121 


151 


158 


151 


N obs 


2575 


305 


73 


27 


fJ'SM 


2436 


249 


56.4 


18.8 


F 


35% 


6% 


41% 


56% 



Table 2: Results of the significance analysis on the M e j spectra: P m i n is the 
minimum probability along the M e j spectra, as defined in Eq. (QJ); N b s is the num- 
ber of observed events at the corresponding M e f, fi$M is the expectation from the 
SM background; E is the fraction of the simulated experiments with the minimum 
probability P < P m ,i n . 





No cos 9* cut 


cos#* <0.4 




e + p 


e p 


e + p 


e p 


p . 

1 mm 


1.8% 


9.9% 


1.1% 


24.5% 


M uj (GeV) 


246 


225 


269 


217 


Nobs 


10 


22 


2 


9 




4.5 


16.2 


0.16 


6.8 


F 


76% 


77% 


24% 


94% 



Table 3: Results of the significance analysis on the M V j spectra. For details, see 
the caption to Tabled 
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Resonance 
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e + u 


V e + U 
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e + u 


S e +d 
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e + d 
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2/3 


e + d 






vu 






vu 


<7 

u e u 


-1/3 


e"u 


Ve-u 


-1/3 


e~u 






vd 






vd 


S e -d 


-4/3 


e~d 


Ve-d 


-4/3 


e~d 



Table 4: First- generation scalar and vector resonant states that can be produced 
in e ± p scattering. The top half of the table lists color-triplet states with fermion 
number F = L + 3B = 0, while the bottom half lists those with F = 2. The left and 
right sets of columns list scalars and vectors, respectively. 
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Table 5: The F = (upper part) and F = 2 (lower part) leptoquark species of 
the Buchmuller-Ruckl- Wyler model ^ and the corresponding couplings. Those LQs 
that couple only to neutrino and quark and therefore could not be produced at HERA 
are not listed. The LQ species are classified according to their spin (S for scalar 
and V for vector), their chirality (L or R) and their weak isospin (0, 1/2, 1 ). The 
leptoquarks S and V differ by two units of hypercharge from S and V, respectively. 
In addition, the electric charge q of the leptoquarks, the production channel, as well 
as their allowed decay channels assuming lepton-flavor conservation are displayed. 
The nomenclature follows the Aachen convention JflOtf . 
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LQ Type (F=0) 




yR 


yR 


V? 


oL 
D l/2 


cR 
S l/2 


dL 
°l/2 


M(GeV)(A=0.1) 


266 


268 


282 


290 


282 


282 


269 


M(GeV)(A=0.3) 


386 


287 


305 


367 


308 


303 


286 


LQ Type (F=2) 


oL 

°o 


Sq 


s R 


5f 


V l/2 


yR 


V L 


M(GeV)(A=0.1) 


276 


273 


248 


275 


248 


274 


273 


M(GeV)(A=0.3) 


351 


298 


273 


300 


277 


302 


313 



Table 6: Mass limit of the U BRW LQs at \=0.1 and 0.3. 
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(c) (d) 



Figure 1: Diagrams for ep scattering at HERA via (a) photon and (b) Z° exchange 
(NC) and via W exchange (CC). The leptoquark diagrams for the same initial and 
final states are c) s-channel LQ production and d) u-channel LQ exchange. Here 
e stands for both e + (positron) and e~ (electron), and v for both v (neutrino) and 
v (antineutrino). 
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Figure 2: Comparison of distributions from 94-00 e + p (a-d) and from 98-99 
e~p data sets (e-h) with the corresponding NC MC samples: (a,e) total transverse 
energy, E T ; (b,f) E — P z ; (c,g) electron tranverse momentum, P^, and (d,h) jet 
transverse momentum, P|>. 
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Figure 3: (a) Comparison of the e + p samples (dots) and the NC SM expectations 
(solid histogram) for the reconstructed invariant mass M e j in the e + p — > e + X 
topology. The data (open squares) and the SM expectations (dashed histogram) 
after the cos 8* <0.4 cut are also shown, (b) The ratio between the data and the SM 
expectation before the cos 6* cut. The shaded area shows the overall uncertainties 
of the SM MC expectation. 
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Figure 4: (a) Comparison of the observed e~p samples (dots) and the NC 
SM expectations (solid histogram) for the reconstructed invariant mass M e j in the 
e~p — > e~X topology. The data (open squares) and the SM expectations (dashed 
histogram) after the cos 8* <0.4 cut are also shown, (b) The ratio between the data 
and the SM expectation before the cos 6* cut. The shaded area shows the overall 
uncertainties on the SM MC expectation. 
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Figure 5: Comparison of distributions from 94-00 e + p (a-d) and from 98-99 e~p 
data sets (e-h) with the corresponding CC MC samples for the selected distributions: 
(a,e) missing transverse energy, P T ; (b,f) E — P z ; (c,g) neutrino energy, E u , and 
(d,h) jet transverse momentum, P^. 
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Figure 6: (a) Comparison of the observed e + p samples (dots) and the CC SM 
expectations (solid histogram) for the reconstructed invariant mass M U j in the 
e + p — > vX topology. The data (open squares) and the SM expectations (dashed 
histogram) after the cos9*<0.4 cut are also shown, (b) The ratio between the data 
and the SM expectation before the cos 6* cut. The shaded area shows the overall 
uncertainties on the SM MC expectation. 
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Figure 7: (a) Comparison of the observed e~p samples (dots) and the CC SM 
expectations (solid histogram) for the reconstructed invariant mass M U j in the 
e~p — > vX topology. The data (open squares) and the SM expectations (dashed 
histogram) after the cos8*<0.4 cut are also shown, (b) The ratio between the data 
and the SM expectation before the cos 8* cut. The shaded area shows the overall 
uncertainties on the SM MC expectation. 
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Figure 8: The constant-X limits as a function of the branching ratios into eq and vq 
(shown on the left and right axes, respectively) , and of the resonance mass (x axes), 
for the scalar resonant states listed in Table^ The dotted line corresponds to limits 
set with only the ep —>■ vX data set, the dashed line with only the ep —>■ eX data 
set; the solid black line is the limit set using both data sets, assuming (5 eq + j3 uq = 1. 
For each limit curve, the area to the left of the curve is the excluded region. Results 
for the resonant states are shown for constant limit of X = 0.1 and X = 0.3. The 
shaded area in each plot shows the area excluded by the D0 experiment. 
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Figure 9: The constant-X limit contour for the vector resonant states listed in 
Table^ For details, see the caption to Fig. [$1 
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Figure 10: (a,b) Limits on the total production cross section for a narrow 
scalar resonant state and (c,d) the corresponding limits for a narrow vector resonant 
state as shown in Table^ Limits derived from NC (CC) data samples assume a 
branching ratio f3 eq ((3 uq ) of 1/2, while the combined eq + vq limits assume branching 
ratios of f3 eq = (3 vq = 1/2. 
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Figure 12: Coupling limits as a function of LQ mass for scalar (a) and vector 
(b) F=0 BRW LQs and for scalar (c) and vector (d) F=2 BRW LQs. The areas 
above the curves are excluded. 
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